Defect-free surface states in modulated photonic lattices 
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We predict that interfaces of modulated photonic lattices can support a novel type of generic 
surface states. Such linear surface states appear in truncated but otherwise perfect (defect-free) 
lattices as a direct consequence of the periodic modulation of the lattice potential, without any 
embedded or nonlinearity-induced defects. This is in a sharp contrast to all previous studies, where 
surface states in linear or nonlinear lattices, such as Tamm or Shockley type surface states, are 
always associated with the presence of a certain type of structural or induced surface defect. 
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Interfaces separating different physical media can sup- 
port a special class of transversally localized waves known 
as surface waves. Linear surface waves have been studied 
extensively in many branches of physics [H. For exam- 
ple, electro-magnetic waves localized at the boundaries 
of periodic photonic structures, such as waveguide ar- 
rays or photonic crystals, have been extensively analyzed 
theoretically and experimentally. The appearance of lo- 
calized surface waves in photonic structures is commonly 
explained as the manifestation of Tamm or Shockley type 
localization mechanisms 

Hi a, 

being associated with 
the presence of a certain type of surface defect. Tamm 
states were first identified as localized electronic states at 
the edge of a truncated periodic potential [3], and then 
they were found in other systems, e.g. at an interface 
separating periodic and homogeneous dielectric optical 
media [5, 6]. 

In discrete systems, such as arrays of weakly coupled 
optical waveguides [7] , different types of linear and non- 
linear states localized at and near the surface have also 
been analyzed extensively. It was found that Tamm sur- 
face waves can exists at the edge of an array of opti- 
cal waveguides when the effective refractive index of the 
boundary waveguide is modified above a certain thresh- 
old |8|, [9|, [lOl, lUl, ll2|, ll3|, [l^ , whereas surface localization 
was considered to be impossible when all waveguides are 
exactly identical, as sketched in Fig. [TJa). In the lat- 
ter case, the beam launched into array delocalizes due 
to diffraction [Fig. [T](b)], and it is also strongly reflected 
from the boundary as illustrated in Fig. [Tfc). 

In this Letter we predict, for the first time to our knowl- 
edge and contrary to the accepted notion, that novel 
type of generic defect-free surface waves can exist at the 
boundary of a periodic array of identical optical waveg- 
uides, which axes are periodically curved along the prop- 
agation direction as schematically shown in Fig. [T][d). 
The periodic bending of waveguide axes was shown to 
result in the modification of diffraction 



for particular values of the bending amplitude, and this 
effect is known as dynamic localization or beam self- 
collimation [l5|, 0, [S]- Under such very special 
conditions, the beam experiences periodic self-imaging, 
propagating without spreading for hundreds of free-space 
diffraction lengths, as illustrated in Fig. [D^e). On the 
other hand, if the beam is launched at the edge of a semi- 
infinite modulated lattice tuned to the self-collimation, 
one can intuitively expect that it can not penetrate deep 
into the lattice as away from the lattice edge effect of 
the boundary is negligible and coupling between the lat- 
tice sites is canceled in the self-collimating lattice. In 
Fig. [iff) we indeed observe that the beam remains lo- 
calized at the surface of the self-collimating modulated 
lattice. However, our most nontrivial finding detailed be- 
low is that surface localization is possible for an extended 
range of structural parameters even when diffraction is 
non- vanishing. 

We study propagation and localization of light in a 
semi-infinite one-dimensional array of coupled optical 
waveguides, where the waveguide axes are periodically 
curved in the propagation direction z with the period 
L, as shown schematically in Fig. [TJd). When the tilt 
of beams and waveguides at the input facet is less than 
the Bragg angle, the beam propagation is primarily char- 
acterized by coupling between the fundamental modes of 
the individual waveguides, and it can be described by the 
tight-binding equations taking into account the periodic 
waveguide bending lH, [l^ , 
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which strength nontrivially depends on the waveguies 
bending and optical wavelength. An interesting fea- 
ture is that the diffraction can be completely suppressed 



Cexp [-ixo(2;)] an+i + Cexp [i±o{z)] an-i = 0, 

(1) 

where an(^) is the field amplitude in the n-th waveguide, 
n = 1, . . ., and an<o = due to the structure termination. 
Transverse shift xo(^) = xo(^ + L) defines the periodic 
longitudinal lattice modulation. Coefficient C defines 
the coupling strength between the neighboring waveg- 
uides, it characterizes diffraction in a straight waveguide 
array with xq = [2o| [see an example in Fig. [IJb)]. 
Expression ([T]) shows that the effect of periodic lattice 
modulation appears through the modifications of phases 
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FIG. 1: (a-c) Beam propagation in a straight semi- infinite lat- 
tice shown schematically in (a): (b) discrete diffraction away 
from the lattice boundary and (c) diffraction and reflection 
from the surface when beam is coupled to a boundary waveg- 
uide, (d-e) Beam dynamics in a sinusoidally modulated lat- 
tice, shown schematically in (d). The modulation amplitude 
is chosen to obtain self-collimation regime (A = Ao — 1.24 for 
the period L = 3.25). Diffraction is suppressed and beam re- 
mains localized at the input location either (e) away or (f) at 
the lattice boundary. 



of the coupling coefficients along the propagation direc- 
tion z. In order to specially distinguish the effects due 
to diffraction management, we consider the light propa- 
gation in the waveguide arrays with symmetric bending 
profiles, since asymmetry may introduce other effects due 
to the modification of refraction, such as beam dragging 
and steering [2l|, [H, [2^. Specifically, we require that 
xo(^) = i{z — z) for a given coordinate shift where 
function i(z) is symmetric, i{z) = f{—z). 

In order to analyze light propagation near the sur- 
face of a semi-infinite modulated lattice, we first con- 
sider the case of small modulation periods L, such that 
the parameter = 27r/L is large, 1. Then we can 
employ the asymptotic expansion (see, e.g., Ref. [HI) 
a.n{z) = Un{z) + ^^^oyn,miz) exp{imKz), where Un(2;) 
have the meaning of the averaged field values over the 
modulation period, and we take into account first- and 
second-order terms for the oscillatory corrections which 
have zero average, Vn,m = Vn,m^^^/^~^ + Vn,m^^^/^~^- Since 
the modulation is periodic, we can perform Fourier ex- 
pansion of the coupling coefficients as C exp[— ixo(^)] = 
Cm exp(im/^z). Then, in the regime close to self- 
collimation when the average coupling is small, |Col ~ 
0{hi~^)^ we combine the terms of the same orders [2^ 
and finally obtain the effective equations for the slowly 



varying functions u.n{z)^ 

. dUn ^ CoUn+1 + CoUn-1 

dz 

+ (5l,n (AiUi + A2U2) + (^2,nA2Ui = 0. 

Here 5 is the Kronecker delta. 
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the bar stands for the complex conjugation, and Un<o = 
0. From these equations one can see that the effect of pe- 
riodic modulation is to introduce the "virtual" defects Ai 
and A2 at the lattice boundary. We now seek solutions in 
the form of stationary modes, Un{z) = Un(0) exp(ikz/L), 
where k is the Bloch wave- number. The values of |k| < 
2 1 Co I L correspond to a transmission band^ where the 
modes are infinitely extended. On the other hand, the 
modes can become localized at the surface of semi-infinite 
modulated lattices is |k| > 2|Co|L, and we find that such 
solutions exist if the modulation parameters are suffi- 
ciently close to the self-collimation condition where |Co| 
is small. Specifically, there exists one surface state if 
a2 — ai < 2 < Qfi + 0^2, and two surface states emerge if 
0^2 — 0^1 > 2, where ai = |Ai/Co| and 0^2 = |1 + A2/C0P. 

We note that if the modulation is symmetric, such that 
xo(2: + L/2) = — xo(2:), then |Cm| = |C_m| and accord- 
ingly Ai = 0, meaning that the modes should always ap- 
pear in pairs. Moreover, this conclusion is valid even be- 
yond the applicability of the asymptotic expansion, since 
we identify the exact symmetry of the model Eq. ([T|) in 
case of symmetric modulations: for each solution an(^), 
bn(^) = (— l)"an(^ + L/2) is also a solution. Therefore, 
in symmetric structures surface modes always appear in 
pairs with the Bloch wavenumbers of the opposite sign, 

k = ±L\Co\{d + d-i), where d = (|1 + A2/C0P - l)'^', 
As an example, we further consider a sinusoidal mod- 
ulation function of the form xo(2;) = A [cos (27rz/L) — 1], 
similar to the one which has recently been employed to 
demonstrate dynamical localization in modulated waveg- 



uide arrays p^, Il9[. In this case, the Fourier coeffi- 
cients can be calculated analytically. Cm = CJm[^A/Ao], 
and Jm is the Bessel function of the first kind of order 
m. The modulation amplitude Aq corresponds to the 
self-collimation condition [HI, [25[, Aq = ^L/(27r), where 
^ 2.405 is the first root of the Bessel function Jq. Since 
the sinusoidal modulation is symmetric, then for each 
modulation amplitude A such that A~.^ < A < A^-^, 
where A~-^ and A^-^ are the left and right mode cut- 
offs, respectively, there exists (at least) a pair of surface 
modes. We use our asymptotic analysis to estimate the 
cut-off values in the case of small modulation periods, 
Af.JAo ^ l±2A2C2Ji[e]J2[e]/(?'(\/2Tl)). 
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FIG. 2: (Color online) (a) Defect-free surface modes in sinu- 
soidally modulated lattice with modulation period L = 3.25. 
Circles and solid lines show modes Bloch wave numbers k 
calculated numerically and using asymptotic expansion, re- 
spectively. Shading marks transmission band of the lattice, 
(b-c) Numerically calculated modes profiles at the input, and 
(d-e) their propagation dynamics are shown for the two com- 
plementary modes marked 1 and 2 in (a), respectively. 



In order to confirm our analytical results, we calculate 
numerically the mode spectrum of the original Eq. ([1]). 
In Fig. [2](a) one can see that for sufficiently small mod- 
ulation periods L there indeed exists a pair of symmet- 
ric surface modes outside the lattice transmission band, 
and the wave numbers of surface modes calculated us- 
ing asymptotic expansion are in excellent agreement with 
those calculated numerically. At the cross-section z = 0, 
one mode has unstaggered input profile [Fig.[2]^b)], while 
the other one exhibits staggered structure [Fig. [21(c)]. We 
note that there is very weak additional phase modulation, 
Im[an]/Re[an] ~ 10~^, in agreement with the asymptotic 
analysis predicting real profiles up to second-order cor- 
rections. In all the figures, we put C = 1, since results can 
be mapped to the other coupling values using a simple 
transformation an (2;, C, L, A) = an(Cz, C = 1, CL, CA). 

We further demonstrate that defect-free surface modes 
in modulated lattices can be effectively generated using 
single-site excitation of the edge lattice waveguide, if the 
lattice modulation amplitude A is between the left and 
the right cut-offs A~ and A^-^ . An example of such sur- 
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FIG. 3: Beam propagation in two modulated semi-infinite lat- 
tices with two different modulation amplitudes and the same 
modulation period L = 9.75 . Top: Lattice modulation is 
larger than the left surface modes cut-off, A ^ 0.934Ao. (a) 
Beam diffracts when launched away from the surface, (b) Sur- 
face wave is generated when the edge waveguide is excited. 
Bottom: Modulation amplitude is less than the left surface 
modes cut-off, A O.8A0. Beam always diffracts, whether it 
is launched (c) inside the lattice or (d) at the surface. 



face wave excitation after some initial radiation is shown 
in Fig. [3l[b), where even though the lattice modulation 
is very close to the left cut-off, A 1.0065A~-^, the sur- 
face wave is still very well localized. In contrast, when 
the beam is launched far away from the surface, it always 
diffracts if A 7^ Aq, as shown in Fig. [31(a). This illustrates 
the fundamental difference between the dynamical local- 
ization in infinite modulated lattices [15, 25], and forma- 
tion of the defect-free surface modes in truncated mod- 
ulated lattices. While dynamical localization is a purely 
resonant effect which takes place just for one single value 
of the modulation amplitude A = Aq [see Fig. lU^e)], the 
families of defect-free surface modes always exist in a fi- 
nite range of the modulation amplitudes sufficiently close 
to the self-collimation value Aq. If the deviation of the 
modulation amplitude from the self-collimation value is 
greater than the one determined by the left and the right 
cut-offs, the defect free modes disappear, and the beam 
always diffracts irrespectively of its input position in the 
semi-infinite modulated lattice, see Figs. [31(c) and[3l(d). 

For large modulations periods the asymptotic analy- 
sis is not valid, and we use numerical simulations to find 
families of the defect-free surface modes. These results 
are summarized in Fig. [41(a), where hatched is the do- 
main of the existence of the defect-free surface modes 
on the (L, A) parameter plane. For small modulation 
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FIG. 4: (Color online) (a) Solid lines bound domain [hatched] 
of existence of defect-free surface modes in the sinusoidal 
modulated lattice calculated numerically. Dashed lines show 
modes cut-offs obtained from the asymptotic expansion. Solid 
shading marks the region where localized modes can not exist, 
(b-c) Families of defect-free surface modes for the modulation 
periods L = 9.75 and L = 19.5, respectively. Shadings mark 
lattice transmission bands. 



periods, the asymptotic expansion provides an estimate 
for the surface modes cut-offs (dashed lines). When 
the modulation period grows, the number of defect-free 
modes increases, as shown in Figs. HJ^b) and (c). For 
the large modulation periods the domain of the exis- 
tence of the defect-free surface modes is basically lim- 
ited by the region where lattice transmission band ex- 
tends to the whole Brillouin zone from — tt to — tt, and 
therefore localized states cannot exist. The region where 
localized modes can not exist [shown with solid shading 
in Fig. ID^a)] is given by the relation L > 7r/(2|Co|) = 
7r/(CJ„,[CA/AJ). 

We note that although the defect-free surface states 
were introduced here for modulated photonic lattices, 
such novel type of surface modes may also appear in 
other fields where wave dynamics is governed by cou- 
pled Schrodinger type equations ([T]) with z standing for 
time. In particular, by introducing special periodic shift 
of lattice potential it may be possible to observe pecu- 
liar surface localization in Bose-Einstein condensates. On 
the other hand, our results indicate the possibility for 
novel mechanism of surface localization of charged par- 
ticles in complex time- varying driving electric fields, for 
which the possibility of the dynamical localization has 
been suggested earlier (iBj . 

In conclusion, we have demonstrated, for the first time 



to our knowledge, that interfaces of modulated photonic 
lattices can support a novel type of generic defect-free 
surface states. Such surface states appear in truncated 
but otherwise perfect (defect-free) lattices as a direct con- 
sequence of the periodic modulation of the lattice poten- 
tial, without any embedded or nonlinearity-induced de- 
fects. This is in a sharp contrast to all previous studies, 
where surface states in linear or nonlinear lattices, such 
as Tamm or Shockley type surface states, are always as- 
sociated with the presence of a certain type of surface 
defect. Using both asymptotic expansion technique and 
numerical simulations, we presented detailed analysis of 
the different families of the defect-free surface states in 
modulated lattices. 
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